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INTRODUCTION 

A  wide  selection  of  literature  is  devoted  to  the  problem  of  diffraction 

of  radio  waves  around  the  surface  of  the  Earth;  a  review  of  the  most  recent 

(2) 

studies  can  be  found  in  Academician  B.A.  Vvedenskij’s  article  .  Interest 
in  this  question  is  justified  in  that,  at  short  distances  of  the  order  of  a 
few  hundred  kilometres,  refraction  of  radio  waves  in  the  ionized  layers  of 
the  atmosphere  can  be  ignored,  while  diffraction  pl3ys  the  decisive  role  in 
radio-wave  propagation. 

Despite  the  fact  that  a  firm  solution  of  the  question  of  diffraction 
from  the  globe  has  already  been  known  for  some  decades,  a  practically 
applicable  approximate  solution  has  not,  up  to  now,  been  obtained.  In  the 
present  work  we  intend  to  fill  this  gap. 


I. 


STATEMENT  OK  THE  PROBLEM  AND  ITS  SOLUTION  IN  SERIES  -FORM 


By  r,  0  and  if  we'  mean  spherical  coordinates  with  their  origin  at  the 
centre  of  the  Earth.  If  we  ignore  the  irregularities  of  the  Earth’s 
surface,  then  its  equation  will  take  the  form  r  =  a,  where  a  =  Earth's 
radius.  Let  us  place  a  vertical  electric  dipole  (  b  >  a)  at  the  point 
r  =  b,  9  =  0.  Rejecting  in  the  field  components  the  time  dependent  factoi 
e~lwt,  we  can  express  it  by  the  Hertz  function  U,  which  will  depend  only 
r  and  on  9.  Having  by  -k  denoted  the  absolute  value  of  the  wave  vector,  y ' 
shall  have,  for  the  field  in  air,  the  expression 


(1,01) 


while  the  remaining  components  of  the  field  will  equal  zero.  Similar 
expressions  will  apply  for  the  field  within  the  earth.  The  U  function 
must,  where  r  >  a,  satisfy  the  equation 


AU  +  k?tJ  =  0 


(1,02) 


and  the  condition  of  radiation  to  infinity 
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r  -v  ro 
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=  0 


(1,03) 


If  b  s  a,  in  sucii  a  way  that  the  source  (dipole)  is  above  the  Earth's 
surface  and  not  on  the  surface  itself,  then  at  Point  r  =  b ,  9  =  0  the  U 
function  must  have  the  property: 
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ikR 
e _ 

R 


+  U*, 


(1,04) 


where 


R  =  ,/  r2  +  b2  -  2rb  cos  § 


(1,05) 


is  the  distance  from  the  source,  while  the  U*  function  remains  finite 

when  kR  -►  0.  On  the  Earth's  surface  the  Hertz  function  U  must  satisfy 

the  limiting  conditions  guaranteeing  continuity  of  the  components  E  and 

H.  at  the  interface. 

<P 

If  we  denote  the  Hertz  function  inside  the  Earth  as  U2  then  these 
limiting  conditions  will  take  the  form 


k2U  =  k22U,>>  jrjr  (r0)  =  (ru2)  where  r  =  a  .  (1,06) 


The  Uo  function  must,  when  0  <_  r  £  a  (in3ide  the  Earth),  satisfy  an 
equation  similar  to  (1,02)  and  remain  finite. 

The  value  k2,  introduced  in  (1,06)  and  later  formulae,  is  determined 
from  the  equation 


k22  =  ek2  +  i~ 


(1,07) 


and  the  condition  lm(k2)  >0.  In  place  of  the  conductance  of  the  Earth  o 
it  is  convenient  to  bring  in  length  l ,  characterizing  the  Earth's  specific 
resistance,  having  assumed 


c 

4no 


(1,08) 


For  sea  water  the  value  of  l  varies  from  0.05  cm  (very  salt  water) 
to  0.5  cm  (slightly  salt  wa‘  »-r) ;  for  the  earth  this  length  is  hundreds 


and  thousands  of  times  greater..  If  we  bring  in  the  complex  dialectric 
constant  of  the  Earth 


then 


(1,09) 


lc2  =  k  /~n  . 


(1,10) 


A  solution  of  the  problem  posed  in  the  form  of  series  is  well  known; 

Here  we  inscribe  the  necessary  formulae  without  dwelling  on  their  derivation. 

We  assume 


*n(x)  = 


!X 
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Wx>1 


<n(x) 


»  / 


HJL 
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<Vx)' (1,u) 


where  JV(y)  is  the  Bessel  function,  while  H^(x)  is  a  first  order  Hankel 
function.  These  functions  are  connected  by  the  relationship 


'J'n  (x)r»^(x)  "  #n<x)5n(x)  ~  (1,12) 


Additionally,  we  introduce  the  particular  specification  for  the 
logarithmic  derivative  from  iji  (x)  > 


x„(x)  = 


*n<x> 


(1,13) 


As  is  seen  from  (1,01),  the  field  at  the  Earth's  surface  is 
expressed  by  the  values: 


U  =  (U; 


r-a 


U  * 
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Ar 


(no 


r=a 


(1,14) 
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For  these  values  we  obtain  the  series: 


(1,15) 

a,  i6) 

ordered  according  to  LeGendre's  polynomials.  Our  task  consists  of  an 
approximate  summation  of  these  series. 

2.  SUMMATION  FORMULA 

The  sums  which  we  have  to  calculate  have  the  form 


1 

kab 


(2n  +  lk(kb) 
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n=o  ^’(ka)- ^n(k2a)?n(ka) 


P  (cos  b); 
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(2n+l) t  (kb) y  (k  a) 

rr  I  - nr — — - p  <«>«  8), 

2  n-o  c  '  (ka)-  j^xn(k2a)Cn(kcl)  n 


S  =  £  uq>(v)  P  i  (cos  b).  (2  ,01) 

1  3- 


In  the  Equ,(l,14),  the  function  <p(v)  with  an  accjracy  approaching  the 
factor  equal  to 


Wkb) 

q>(v)  *  - j-*- — -  (2,02) 

?v-J  ~  ^7  X'  "i  ^k2G^v-J  ^ka^ 

In  Equ.  (1,15)  it  differs  from  (2,02)  by  the  factor  xv_j(k3a). 

For  a  direct  calculation  of  the  sum  it  would  be  necesssary  to  take 
the  number  of  terms  approximately  equal  to  2ka,  i.e.  equal  to  double  the 
number  of  waves  present  on  the  Earth's  surface.  Since  this  number  is 
enormous,  it  is  evident  that  an  immediate  summation  is  impossible. 

To  compute  the  sum  S  it  is  necessary,  having  taken  advantage  of  the  fact 
that  <p(v)  is  an  analytical  function,  to  convert  this  sum  to  an  integral  which 
can  be  computed  by  one  or  other  of  the  approximation  methods.  A  similar 
conversion  was  first  proposed  by  Watson  in  1918  and  adopted  by  many  authors. 
However,  all  the  authors  attempted  to  apply  the  converted  expression  to  the 
sum  of  deductions,  wl areas  our  intention  is  to  separate  out  the  main  terms 
most  easily  subject  to  investigation  and  to  estimate  the  remainder.  In  this 


method  Oi  omputation  the  main  terms  are  not  predetermined. 

For  the  completion  of  our  conversion  it  was  essential  to  keep  in 
mind  the  following  general  properties  of  the  function  <p(v).  It  is  an 
analytic  function  from  v,  meromorphic  in  the  right  half-plane;  it  has 
poles  only  in  the  first  quarter,  while  in  the  fourth  quarter  it  is 
holomorphic.  At  infinity  ip(v)  declines  so  rapidly  that  all  -the  integrals 
being  investigated  converge. 

As  seen  in  Equation  (2,01)  the  LeGendre's  functions  can  be  expressed 
by  the  function 
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2  sin  d 


(2,03) 


where  F  is  the  sign  for  the  hypergeometric  function, 
denote  the  result  of  a  change  in  G  and  in  P  .  =  P  , 

V  V“j  V—J 

value  0  in  '(  -  ?.  We  shall  then  have 


By  G*  and  P*... 
v  V“J 

(cos  0)  of  the 
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v-J 


_ !_ _  ["ei"*‘l7  0*  +  e"lv8+ir  Cv] 

7i /  2  sin  0 


(2,04) 


As  can  be  seen  from  (.2 ,03),  if  v  lies  within  a  certain  sector  includ¬ 
ing  the  negative  axis,  and  if  Iv  sin  $|  is  large,  the  function  G^  (and 
also  G*)  approximately  equal 


G 

v 


(2,05) 


Substituting  (2,05)  in  (2,04),  we  obtain  the  known  asymptot 
sion  for  P^_j.  if  B(v)  denoted  the  first  term  in  Equation  (2,04) 
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(2,06) 


one  can  substantiate  the  equation 


i(v—  })n 

P*  .  =  e  P  .  +  2i  cos  vtiB(v) 

v-J  v-J 


(2,07) 


which  we  shall  make  use  of  later.  We  note  that  the  function  B(v)  is 
holomorphic  in  the  right  half-plans. 

We  examine  in  the  plane  of  the  complex  variable  three  contours. 
Firstly,  the  loop  Cj  around  the  origin  of  the  coordinates  alleviating 
the  positive  material  axis  and  essential  in  the  positive  direction 
(anti-clockwise).  Secondly,  the  broken  line  C2  encompassing  the  first 
quadrant  ant  running  from  left  to  right  (in  its  horizontal  part  passing 
a  little  higher  than  the  material  axis).  Thirdly,  the  straight  line  C3, 
passing,  through  the  origin  of  the  coordinates,  lying  in  the  second  and 
fourth  quadrants  and  inclined  at  a  small  angle  to  the  imaginary  axis; 
this  straight  line  runs  from  top  to  bottom. 

We  can  write  the  sum  S  in  the  form 
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1 
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j  vip(v) 
Cl 


sec  vtiP*  ,dv, 
v-J 


(2,09) 


since  the  integral  on  the  right  comes  to  deductions  at  the  points 
’  =  n  +  y  .  Since  the  function  <?(v)  is  holomorphic  in  the  fourth  quadrant, 
we  can  substitute  Contour  Cj  for  Contours  C2  and  C3  and  write 


S  -  -  ±  J  vf(v)  sec  vnP*v_jdv  +  ^  j  v(p(v)  sec  vnP*v_^dv.  (2,09) 

C2  C3 

The  normal  conversion  of  the  sum  comes  down  to  this:  the  integral 

at  C3  is  assumed  to  equal  zero  because  of  the  smallness  of  the  uneven 

part  of  (v )  (its  definition  will  be  given  below), and  the  integral  at 

C,>  boils  down  to  the  sum  of  deductions,  but  we  shall  take  a  step  further 

and  break  up  C,.  into  main  term  and  correction  [term].  Substituting  in 

this  integral  the  expression  (2,07)  for  P*  , ,  we  shall  have 

v  “} 


S  =  Si  +  So  +  s3, 


(2,10) 


Sj  =  vip(v)B(v)  dv, 


(2,  IT) 


-if 


/  \  iv71  ,»  i 

v<p(v)  sec  vnc  P  dv, 

v—5 


(2,12) 


S3  - 


v<p(v)  sec  vnP*  _jdv. 


(2,13) 


In  Integral  S(  the  subintegral  function  is  alieady  without  poles  on 
the  material  axis  (and  also  in  the  fourth  quarter),  therefore  Contours 
C?  and  Cp  [sic]  for  it  are  equivalent.  We  denote  by  C  any  contour  equivalent 
to  Co  and  C  ] , 

The  presentation  of  the  value  S  in  the  form  of  three  integrals  (2,10)  is 
accurate;  in  the  derivation  nothing  was  neglected.  But  the  definition  of 
integral  Sp  and  S<  indicates  that  these  integrals  are  minute  in  comparison 
with  S 2 0  In  truth,  if  we  are  to  compute  the  integral  as  a  sum  of  deductions 
in  the  poles  <p(v)  then  we  are  convinced  that,  in  relation  to  Sj,  it  will  be 
of  the  order  of 


i’ivj  (n”9) 


(2,14) 


Where  vj  is  closest  to  the  material  axis  of  the  pole  <p(v),  The  imaginary 
part  of  V]  is  positive  and  at  high  values  of  ka  will  be 


Lm(\'i) 


cfka) 


(2,13) 


where  c  :s  a  number  of  the  order  of  unitv  (far  the  absolute  conductor 
c  *  0.70).  Since  ka  is  a  very  large  numlur,  of  the  order  of  millions 


(for  A  =  40  in,  ka  =  1.000,000),  it  is  apparent  that  the  value  (2,15)  will 
he  large  (for  instance  equal  to  70),  whi  le  the  value  (2,14)  will  be  very 
small, (In  our  problem  0  may  not  be  close  to  n ,  since, then,  in  consequence 
of  the  necessity  to  consider  the  effects  of  the  ionized  layers  of  the 
atmosphere,  our  equations  would  generally  cease  to  be  applicable).  As 
far  ns  the  integral  S3  is  concerned,  its  value  is  determined  by  the 
uneven  part  of  the  function  <p(v).  But  the  uneven  part  o£  q>(v)  will  be 
of  the  order  of 


'ik?fl 


(2,16) 


And  since  the  i  .'ginary  part  k2a  is  positive  and  extremely  large, 
then  tne  value  of  *5,16)  will  be  an  inexpressibly  small  number. 

The  following  physical  picture  gives  a  graphic  presentation  of  the 
smallness  of  Integrals  S?  and  S3.  Integral  S2  is  the  amplitude  of  a  wave 
travelling  without  a  break  (as  a  result  of  a  single  diffraction)  around 
the  world  once  or  several  times.  Integral  S3  is  the  amplitude  of  a  wave 
penetrating  right  through  the  Earth  with  the  attenuation  normal  to  the 
Earth.  Obviously,  both  these  integrals  are  minute  in  comparison  with  the 
amplitude  of  a  wave  travelling  from  the  source  through  the  air  by  the 
shortest  path. 

Thus,  with  all  accuracy  attainable  by  the  statement  of  the  physical 
problem,  the  sum  S,  determined  by  the  equation  (2,01),  may  be  expressed 
in  t.he  form  of  an  integral  Sj ,  which,  after  substituting  for 
the  equation  (2,06),  takes  the  form 
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v<p(v)e1,;^G*  dv. 

V 


(2,17) 


3.  CALCULATION  OF  THE  HERTZ  FUNCTION  FOR  THE  ILLUMINATED  REGIONS 

If  we  understand  by  <p(v)  the  equation  (2,02),  then  the  connection 

between  sum  S  and  the  value  U„  will  he 
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U 
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Lab 


S. 


(3,01) 


Therefore  our  approximate  expression  for  is  written 
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nkab  /  2  sin  9 


v(p  (v)e 


iv$ 


G*  dv 
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(3,02) 


The  position  of  the  main  integration  section  in  Equation  (3,02)  will 

depend  on  which  point  of  the  integral  is  being  computed.  Generally  speaking, 

the  main  section  will  lie  close  to  v  n  v  ,  where 

o  * 
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kh 
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_ Qb  sin  0 _ 

y  a?-  +  b2  -  2<2b  cos  d 


(3,03) 


The  value  h  is  the  length  of  the  perpendicular  dropped  from  the  centre  of 
the  Earth  to  a  beam  [of  light]  (i.e.  to  a  straight  line  joining  the  source 
to  the  point  of  observation). 

For  an  approximate  computation  of  Integral  it  is  essential  to  find 
for  G*y  and  <p(v)  the  asymptotic  expressions  which  will  be  applicable  in 
the  main  section'.  Since  the  values  v0  and  vq6  are  large  in  comparison  with 
unity,  we  can,  in  accordance  with  (2,05),  assume 


G* 

v 


(3,04) 


For  the  Hankel  functions  entering  into  cptv)  one  may  attempt  to  use 
Debye's  expression 


where 
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(3,05) 


(3,06) 
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.  ' (it  «•  .s  ittiM  are  ns. •<.’  in  4’ontl  i  I  ions  where 
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•  j. 


A-»  1 .1  r  .is  tin*  function 
expression  uscil  I  if  i  t 


i  ,  (k  .</)  in  concerned,  in  the  proximity  of  v 
with  sufficient  nccurncy,  is 


V 


0,08) 


i,.  ,  xpl.tiu  in  wli.il  coml » t i w*  the  disparity  0,07)  is  used,  wo  denote 
t, v  ,  the  angle  between  the  vortical  at  the  point  of  observation  and  the 
«li  left  i<>n  to  t !,e  source  and  introuuce  the  parameter 


(3,09) 


it  is  easy  to  eo  that  for  v  -  v .  ,  m  •  krt,  the  disparity  (3,07)  is 
equi  v.i  lent  I..  the  condition  such  that  p  is  a  large  positive  number.  Such 
value-  -5  i*  are  applicable  to  the  illuminated  region.  Values  of  p  of  the 
ini.  t  ‘I  unity  (positive  and  negative)  are  app  I  ieab’e  to  the  twilight  region; 
ui,,  p  -  t!  it  gives  the  edge  of  the  geottic trie  shadow  (the  horizon  line). 

u  .;v.  an. I  large  values  in  relation  to  absolute  p  values  are  applicable 

l .  *  f  *n  •  i i  uli* 1 1  r t* i\\  oh  •  <•' 


In  this  paragraph  we  examine  the  condition  of  large  positive  values 
,,|  ()  'ill  uiuinat  ed  region);  other  conditions  will  he  examined  in  subsequent 

paras . 

When  p  1 ,  as  we  have  seen,  Debye  expressions  are  applicable  to  the 
U. ink.- 1  functions.  Substituting  them  in  (3,02)  and  using  (3,04)  and  (3,08) 

w.-  obtain  the  integral 
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(3,1.) 


Ip  eondi lions  where 


kit  cos  y  ■s 


(3,12) 


whcr»*  It  b  -  i  is  the  height  of  the  source  above  the  Earth,  Equation  (3,10) 
can  be  computed  according  to  the  fixed  phase  method,  obtaining  the  "reflect¬ 
ing  fo  nmi  In" 
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itcR 
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(3,13) 


who  re 


R  =  \  <r‘  +  b‘‘  -  2«b  cos  0 


(3,14) 


is  the  distance  i rom  the  source,  while  W  is  the  "attenuation  function", 
which  in  our  case  is  equal  to 


W  = 


1  + 


(3,15) 


l  - 


sin' 


sec  y 


The  value  of  p  * ,  determined  by  the  series  (1,16),  differs,  in  our 
approximal ion ,  only  by  the  constant  multiple  from  11  ,  namely 

ii 


b 

•  t 


(3,16) 


l he  last  equation  is  correct,  not  only  lor  the  illuminated  region, 
but  in  all  cases. 
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II  C"<>  -n>n  0,12)  is  not  achieved,  then  the  denominator  and  subin- 
i  gr.il  equation  (3,10)  will  not  be  a  slowly  changing  function.  If  we 
as  unite  that  tiic  conditions  are  fulfilled. 


1  •'  ~7  * '  (ka)  (3,17) 

ir 

I  kR  •  yy  (3,18) 


(the  lesult  ol  which  is  the  disparity  p  > *  I),  thus  the  equation  (3,10) 
(.in  he  computed  appro  si  mat  e  l-y  by  introducing  a  new  variable 
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(3,19) 


for  the  function  W  in  (3,13)  the  approximate  expression 


W  = 


_  .kR,  . 

y  2kR  |  ti“L  2  (l‘  l,  ,) 


l‘d|i . 


(3,20) 


is  obtained,  where 


(3,21) 


is  the  inclination  ol  the  ray  to  the  horizon  and  Contour  F  is  a  straight 
line  p.issinc,  through  the  point  u  »  u,,  from  the  quarter  in  the  second 
quarter  o|  the  plane  e  (more  precisely  n  -  p,,).  The  equation  (3,20)  can  be 
(Oinputed  without  further  neglect  and  gives  the  known  Weyl_van  der  Pol 
lo rum  I  a.  II  we  assume 
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(3,22) 


-I  J- 


then  we  sli.ill  have 


n+..  i 

W  =  2  -  [  e^  d«  .  (3. 23) 
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i  •« 

In  nitlcr  I  o  obi'ii  i  n ,  Iron)  our  oppressions  U  and  U  ’,  the  formula  for 

a  <z 

the  lieltl,  ii  is  nooess.iry  to  di  ffer'-.ftisite  these  expressions  according 

i  >  >i.  lUil  i  no  derivatives  according  to  0  are  easy  to  compute,  since,  in 

iicR 

(3,1  i) ,  -it  is  possible  lo  consider  all  multiples,  apart  from  e  as 
i  oust  .111  L  s . 

4  •  ASYMPTOTIC  KXPKI-SSIPN.S  FOR  Till-  HAN  Kill.  FUNCTION 

baler,  we  shall  have  to  investigate  the  case  where  the  observation 
point  lies  in  the  Lwi light  region. 

This  condition  is  characterized  by  the  values  of  the  p  parameter 
(positive  and  negative)  of  the  order  of  unity.  Since  the  disparity  (3,07) 
ceases  to  he  uselul  for  these  values  of  p  in  the  main  section  of  integra¬ 
tion,  the  Delve  expressions  (3,05)  for  the  llankel  function  become  inappli¬ 
cable  and  must  be  exhanged  tor  others.  New  expressions  for  the  llankel 
l unctions,  appropriate  to  our  task,  may  be  obtained  from  asymptotic  expres¬ 
sions  quoted  in  our  former  work^  and  also  from  formulae  in  Watson’s  well 
known  hook  *  , but  it  is  simpler  to  derive  them  direct. 

Unknown  expressions  give  a  description  of  llankel  functions  via  the 
I  unction  w(l),  determined  by  the  integral 


w(t) 


tz-  -  z 


dz , 


<4  ,01) 


where  the  i out  our  l  runs  from  infinity  to  zero  along  the  straight  arc 
/.  -  )  nul  l  rom  zero  to  infinity  along  the  straight  arc  z  =  0  (along 

the  positive  material  axis).  The  function  w(t)  satisfies  the  differential 
equal  i o:i 


w'(t)  =  tw(t) 


(4,02) 
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i>l  (  he  ini!  i  .1 1  rond i  1  i ims 


w(0) 


w'  (0) 


=  1.0899290710  +  i  0.6292708425, 


=  0.  7945704238  •-  i.  0.4587454481. 


r(4 ,03) 


li  is  tlu>  entire  transcendental  function  which  breaks  down  into  a 
step  series  of  the  form 


w(t)  -  w(0)  I  +  2>3  +  (2-5) (3*6)  +  (2*5*8) (3*6*9) 


-} 


(4,04) 


t«.  t7 

*  w’(0)  t  r  -i-  +  - 


.10 


3-4  (3*6) (4*7)  (3-6-9) (4*  7 • 10) 


+  . . . 


II  we  separate  in  w(i)  the  material  and  imaginary  parts  (for  the 
material  l)  ■>nd  assume 


w(t)  =  u(t)  +  iv(t) 


(4,05) 


then  u(l)  and  v(t)  will  he  two  independent  integrals  of  the  equation  (4.02), 
connected  by  the  relationship 


u * ( t ) v ( t y  -  u(t)v* (t)  =  1 


(4,06) 


Asymptotic  exp  1  ess i 011s  lor  these  functions  lor  large  negative  t 
values  aic  obtained  by  separating  the  material  and  imaginary  parts  in  the 


equal  ions 


w(t)  = 


(4,07) 
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w'(f) 


t;v  t)^;/3i(- 1)3^( 


(•4,08) 


For  |arj;o  positive  L  values  asymptotic  expressions  for  u(t),  v(t)  and 
their  derivatives  have  the  form: 

2  Vi  1  2  3/2 

_  1  j  l)  rrt  -r  —  t 

u(t)  =  t  e  ;  u* (t)  =  t  e  ;  (4,09) 


v(l)  - 


I  -1/" 


2  L 


- 1  t7h 

e  J  ;  v'(t) 


1  / 

2  C  e 


1  •  I  t3/? 


.  (4,10) 


From  i ho  series  (4,04)  it  is  not  difficult  to  extract  the  formula 


to 


.  ii 

*3 


L“ 

2e  6  v(-  t), 


(4,11) 


f  ,2n.  ) 
te  3 


.  n 
l3 


lu(t)  -  iv(t) ) ,  (4,12) 


which  j>ive  an  impression  of  the  behaviour  of  w(t)  in  a  complex  plane. 

We  note  that  the  function  w(t)  is  expressed  via  the  Hankel  function 
ol  the  order  1/1  by  the  lormula 

_  .  2ji_  _1_  ^ 

w(i  )  =  /  je  3  (-  t)2  Ha  j  ^  |  (“  t)  ''  j  (4,13) 

T 

llavinn  learned  the  major  properties  of  w(t),  we  switch  to  the  deri¬ 
vation  of  the  asymptotic  expression  for  the  Hankel  function  H^(p), 
where  e  and  ,•  are  larju  ana  close  to  each  other,  so  that  the  relationship 


v  -  I1 
/  l>/2 


(4,14) 


remains  limited. 


-lb- 


Il»t‘  llankei  I  unction  assumes  the  integral  represc»nL;it  ion 


H*i;(p) 

v 


_L  ‘ 

«  i 


-p  sh  v  +  vv  , 
e  dv, 


(4,15) 


where  the  contour  C  {toes  along  the  straight  line  lm(v)  -  -  n  from  -  ni  -  «* 
to  a  certain  point  v  =  v  (e.g.  v  =*  -  n/^~T~  in),  then  along  the  straight 
line  from  v  =  v  to  v  =  0  and,  finally,  along  the  material  axis  from  zero 
to  infinity.  We  express,  in  accordance  witii  (4,14),  v  through  t  and  intro¬ 
duce  the  variable  integration 


■  Y\'- 


(4,16) 


Considering  t  and  z  to  be  finite  and  p  to  be  large,  we  may  displace  the 
subintegral  function  in  (4,15)  according  to  the  negative  (fractional)  power 
of  ,i.  Since  in  the  main  part  the  transformed  S  contour  coincides  with  the  1’ 
contour,  we  can  write 


ii(1)(i')  =-4-  £ 

v  in  2 


J_  fi.  3 

60  2 


!**  +  . . .  "J  dz, 


(4,17) 


and  calculating  the  integrals  with  the  aid  of  (4,01) 


»<!>(*->  —  (?)  'P{"M  -5i  (  r)  '/3“C)(0  *  .  0.: 

/  n  t  *>  '  * 


lly  virtue  of  the  differential  equation  (4,02),  the  5th  derivative  equals 


w^  ‘  u  )  =■  t  w'  (t)  +  4lw(t). 


(4,19) 
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iih:: I  i  t  lit  i ii )•  this  expression  in  (4*18)  and  going  over  by  Formula 

(,.,!!»  (o  (lie  I  unction  (.  J_  (f> ) ,  we  shall  have 

2 

V*(,°  =  "  ‘  (  2  ]1/f’  f (t)  “  6^  (  2-)  ?/3  [(t  W’(C)  +  4tw(t)]+  •*•}. 

(4,20) 

Differentiating  this  expression  according  to  p  when  v  is  constant 
and  taking  into  account  t he  dependence  of  t  on  p,  we  get  the  following 
formula  for  the  derivative: 


We  shall  be  making  use  of  these  expressions  later. 

r) .  EXI* HESS  ION  OF  THE  HERTZ  FUNCTION  APPLICABLE  IN  THE  TWILIGHT  REGION 

We  transcribe  Expression  (3,02)  for  the  Hertz  function,  replacing  in  it 
the  value  0*  by  the  approximate  value  /  ir/y  and  the  value  sin  0  in  front  of 
the  integral  by  the  approximate  value  0.  We  obtain 


U 


a 


|  (p(v)e*V'>  •/r-vdv 
C 


(5,01) 


by  the  0  contour  we  shall  understand  the  C?  contour  described  in  Para. 2, 
or  sp.ne  equivalent  of  it.—  The- main  part  of  ttiir'integratiorT  will ,  in  our  case 
(i.e.  for  finite  values  of  the  p  parameter),  1 ie  close  to  the  point  where 
v  =  ka. 

The  function  xy_j  (k  ><i)  appearing  in  the  expression  (2,02)  for  q>(v),  we 
can  now  replace,  therelore,  by  the  value  of  the  expression  0,08)  when  v  =  kfl, 
alter  which  we  shall  have 


'p(v) 


(5,02) 
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I'or  tin.*  f'unction  and  its  derivative  we  must  take  expressions  which 

will  Ik-  correct  in  the  proximity  of  v  =  ka,  Such  expressions  were  obtained 
in  tlie  preceding  paragraph.  Leaving  the  main  terms  in  (4,20)  and  (4,21), 
we  shall  have 


r  \m  \  1  /  f> 

Vj<ka)  ■  ~  i(  y)  w(t>*  (5,03) 

f.;_-j-(k«)  »  i(  y  )"l/(w’(t),  (5,04) 


where  t  is  linked  to  v  by  the  relationship 


(5,05) 


The  numerator  in  (5,02)  is  obtained  from  (5,03)  by  substituting  a  for  b 
and  t  for  t',  where 


v  =  kb 


(5,06) 


Conij  a  ring  (5,05)  and  (5,06)  we  get  a  link  between  t  and  t'.  But  the 

rel.ti  ionship  h/a  ,  where  h  =  b  -  a,  is  a  small  value  (we  will  consider  it 

-71 

to  b.*  ol  the  same  order  as  (ka)  ).  Ignoring  this  value  in  comparison 

wiilumi  tju.  JKe-.ace  .ahle_t.ft_W.tjlte  .  .  - - -  --- 


t'  -  t  -  y,  (5,07) 

where 

kii 

y  =  -  (5,08) 

(kf  )U" 
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is  a  value  proportional  to  the  height  of  the  source  above  the  Earth.  The 
value  y  may  be  called  the  accepted  altitude  of  the  source.  Thus,  with 
accuracy  to  terms  of  the  order  of  h/a  or  (k a)  2/\  will  be 


( 


.v 

}  w(t  -  y), 


(5,09) 


where  t.  is  deler  ined  from  (5,05)  (in  the  multiple  before  w(t  -  y)  we 
similarly  substitute  b  for  a). 

Substitutions  (5,03),  (5,04)  and  (5,09)  in  (5,02)  give  the  approximate 
expression  for  <p(v). 

If  we  assume  for  brevity 


<1  »  i 


we  shall  have 


'p(v) 


(  h(L  <1 1/n  w(t  -  y) 

(  2  '  w’(t)  1  qw(t )  * 


(5,10) 


(5,11) 


Recalling  the  formulae  (1,09)  and  (1,10)  we  can  write  the  value  q 
in  the  form 


=  1  l  t) 


1/3 


"  1  *  *2«T 


C  +  l 


2  nl 


(5,12) 


or,  with  the  same  accuracy 


<1  *•  i 


■  (?) 


l/'i 


+  l  +  i 


.  X 


m 


(5.13) 


The  last  expression  is  rather  more  convenient  for  calculation. 
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It  remains  for  us  now  to  substitute  the  value  (p(v)  from  (5,11)  in  the 
lonmila  (5,01)  and  transfer  it  to  die  variable  integration  t.  Making  this 
substitution,  we  can  exchange  the  value  /v~  under  the  integral  with  the 
constant  value  /  ka,  and  also  write  a  in  place  of  b  in  the  multiple  before 
the  integral.  As  a  result  we  obtain  a  formula  which  may  be  written  in  the 
form 


U 


a 


ixt  w(t  -  y) 
w' (t)  -  qw(t) 


dt, 


(5,14) 


where  letter  x  signifies,  for  brevity,  the  amplitude 


X  = 


(5,15) 


which  can  be  called  the  applied  horizontal  distance  from  the  source, 
such  that  y  and  q  have  the  value  (5,08)  and  (5,13).  Countour  C  must  embrace 
all  the  poles  of  the  subintegral  function,  as  we  shall  see,  all  of  them  are 
located  in  the  first  quarter  of  plane  t. 

In  order  graphically  to  present  the  relationship  of  the  horizontal 
and  vertical  scales  in  the  variables  x  and  y,  we  shall  compose  an  expression 
for  tlie  parameter  p,  determined  by  the  formula  (3,09).  From  examination  of 
a  triangle  with  its  angles  at  the  centre  of  the  Earth,  the  source  and  the 
poiiw  of  observation,  it  is  not  difficult  to  extract  the  approximate 
to  ram  la 


Y 


(5,16) 


Thus,  the  equation  lor  the  horizon  line  is  x  =  »  y  .  Later  we  shall 
require  a  link  between  the  distance  U  from  the  source,  considered  along  a 
straight  line,  and  the  horizontal  distance  »  h,  i.e,  (ka )  ^  \  »  y, 
this  link  has  the  lorm: 


kR  =  kaO  + 


(5,17) 
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WllO’'<> 


O 


(5,18) 


6.  STUDY  OF  TUI-  EXPRESSION  FOR  THE  HERTZ  FUNCTION 

Hit*  expression  (5,14)  found  for  the  Hertz  function,  may  conveniently 
lx*  written  in  the  form 


ik.tO 

U<t  =  '  V(x‘  y‘  (6, on 


where 


V(t,  v,  q) 


e^c  y(L~  ~  v) _ 

w' (t)  -  qw(t ) 


dt. 


(6,02) 


Tlu*  value  V,  by  analogy  with  tiie  value  W  introduced  earlier  in  (3,14), 
cun  lx*  called  Liu*  attenuation  factor.  We  shall  establish  the  connection 
between  V  and  W.  Since,  in  the  denominators  of  Formulae  (3,13)  and  (6,01) 
we  can  consider  the  values  ot  R  and  uO  as  being  equal  one  to  another,  we  sha‘1 
obtain,  in  consequence  ol  (5,17): 


W 


Ve 


■l«Y 


(6,03) 


We  now  hive  to  examine  tlu*  expression  (6,02)  for  V.  We  shall  begin  this 

xifitfllLLVa  Ujl'i a* .  ill1 '_cu.se  where  the  value  p  i  s  ..post  t  i  vc .  and_  1  urge  (illumina- 

tcil  region).  Wc  have  already  examined  this  case  from  a  different  approach 
(1‘ura.  1).  However,  since  formula  t6,l)2)  war.  brought  out  by  tis  for  finite  p, 
it  is  an  interesting  conclusion  that  it  is  true  also  for  large  values  of  p. 

When  p  *■ I,  the  main  P  irt  ot  the  integration  will  lie  with  the  larger 
negative  t  (namely  close  to  .  =  p).  Using  Kxpressions  (4,07)  and  (4,08) 

lor  w  and  w1  and  using  the  tixed  phase  method,  we  obtain 
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7 


2 


i 


i  -  -A 

p 


(6.04) 


and  *  as  a  result  of  (o,03), 


2  (6,05) 


P 


I'lu*  last  expression  practically  coincides  with  (3,15).  We  note  that  if 
,  is  o|  the  order  ol  unity  or  if  large,  then  for  the  applicability  of  the 
fixed  phase  method  the  condition  p  » 1  suffices;  if  x  is  small  then  it  is 
necessary  that  y'  *s  2 x •  If  the  latter  condition  is  not  fulfilled,  but 
the  condition 


y 


(6,06) 


is,  then  (he  integral  (f> ,02)  may  he  calculated  by  a  different  method. 

In  the  asymptotic  cxpicssion  for  w(t  -  y)  a  further  simplification  may  be 
made,  alter  which  the  integral  (4,02)  takes  the  form 


i-r  , -  r  i  .  t  +  i  y>  -t 

V  «  e  '*  y]  ~  |  - . —  - . .  dt.  (6,07) 

r  *  "  t  "  iq 


Having  taken  as  a  variable  integration  the  value  /  -  t,  we  arrived  at 
an  integral  in  the  lorm  (3,20)  in  which  j  /  -  t  =  (— ]  and  once  again 
get  the  Weyl  -  van  der  Pol  formula  with  values  of  o  and  t,  cquaT~to* 


.  » 

.  J. — 


(6, OH) 


ami  [ii  n  t  1 1  1 1 1  y  idiih  ide  w;  •  h  (1,22). 
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U'c  now  proceed  to  n  more  interesting  condition,  when  the  value  of  {>  is  a 
positive  or  negative  number  of  the  order  of  unity.  As  we  know,  this  is  the 
twilight,  region,  where  the  diffraction  phenomena  plays  a  major  role. 

If  the  value  of  x  and  y  are  of  the  order  of  unity,  then  the  most 
convenient  method  of  computing  the  integral  (6,02)  is  to  present  it  in  the 
form  ol  the  sum  of  deductions  relating  to  the  poles  of  the  scbintegral 
funol i on . 

It  ts  =  ts(q)  signifies  tile  roots  of  the  equation 


w'(t)  -  qw(t )  =  0 


(6,09) 


then  we  shall  have 


V(x,  y.  <i) 


elxt  s  wts  ~  y 
tg  -  q?  w{ts) 


(6,10) 


The  roots  t s (q )  .ire  the  essence  of  the  function  from  the  complex 
parameter  p.  When  q  =  0  they  revert  to  the  root  tg  -  ts(0)  of  the  derivative 
w'(t),  hut  when  q  -  •”  they  revert  to  Root  ts°  =  ts(°').  The  values  tg  and  t° 
have  Phase  n/l,  so  that 


.  IT 


(6,11) 


We  show  here  a  module  of  the  first  few  roots 


s 

U 

It  °i 

1  s  1 

1 

1.01879 

2.33811 

2 

3.24820 

4.08795 

1 

4.82010 

5.52056 

4 

6.16311 

6.78671 

') 

7.37218 

7.94417 
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At  large  s  values  will  be 


(6,12) 


At  linite  values  of  q;,  one  may  compute  the  roots  using  the  differen¬ 
tial  equation 


dr, 

) 

dq 


(6,13) 


which  is  easily  extracted  from  ,4,03).  Root  t  can  be  determined  either  as 
the  solution  to  (A,  13),  which,  when  q  =  0,  reverts  to  t^,  or  as  the  solu¬ 
tion  which,  when  q  =  ••*,  reverts  to  t  °;  both  solutions  coincide.  From  the 
first  definitio  it  is  easy  to  construct  for  ts  a  series  according  to 
rising  degrees  of  q;  it  will  converge  when  |q|  <  /  ts | .  From  the  second 

del i n i t i on  it  is  possible  to  construct  a  series  according  to  falling 
(negative)  degrees  or  q;  it  will  converge  where  |q|  >  |  /~t~  |.  We  are 
not  reproducing  th*»se  series  here.  We  note  that  the  value  which  at  large 
values  of  |q|  is  close  to  q-’,  is  not  a  root  of  Equation  (6,09). 

In  conditions  where  y*  «  2  |  / tg  |  we  have  the  approximate  equation 


w(ts  -  y)  _ 

"w‘(t  )  =  ch  (y  /~T7)  -  ~^=ZS^  (y  '  t  ),  (6,14) 

s  /  t 


which  permits  evaluation  ol  the  remote  terms  of  series  (6,10).  At  large 
i  values,  such  that  |q|  <.<.  |  /~t~  |,  approximately  t  will  equal  t  (0) 
which  will  equal  From  here  and  from  (6,14)  it  follows  that  the  series 

(6,10)  always  converges.  However,  if  x  is  small  or  y  large,  then  to  compute 
the  sum  ol  the  series  may  require  a  large  number  of  terms. 


in  the  shaded  region,  when  the  value  of  p  is  large  and  negative,  the 
seiies  (6,10)  converges  very  quickly  and  its  sum  approximately  boils  down 
to  the  lirst  term. 

Our  series  (6,10)  corresponds  to  the  Watson  series,  but  has  the  advan¬ 
tage  that  the  terms  ol  our  series  have  simple  expressions. 
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Our  basic  formula  (6,02)  permits  investigation  not  only  of  extreme 
conutLions  (large  positive  p  -  the  illuminated  region,  large  negative  p  - 
the  shaded  region)  but  also  the  intermediate  conditions,  in  particular  the 
twilight  region.  While  for  the  extreme  conditions  our  formula  does  no  more 
than  biarify  earlier  known  formulae  (reflecting  formula  and  Weyl  -  van  der 
Pol  formula  for  the  illuminated  region  and  the  Watson  series  for  the  shaded 
region),  for  the  twilight  region  it  gives  e  sentially  new  results. 

Particular  interest  is  presented  by  the  case  whbre  x  ana  y  are  large, 
while  p  is  finite  (short  waves,  twilight).  This  condition  has  not  been 
investigated  by  anyone  to  date,  and  previously  known  formulae  are  inapplicable 
to  it.  We  shall  derive  here  the  approximate  formulae  which  facilitate  its 
study. 

We  introduce  the  value 


z  =  X  ~  >/~y,  (6,15) 

which  is  the  reduced  distance,  calculated  not  from  the  source  but  from  the 
geometric  shadow  boundary.  In  the  region  of  geometric  shadow  z  >  0,  in  the 
"visible"  region  z  <  0.  We  have 


P 


Lli! 

2X 


z 


+ 


z2 

2X 


(6,16) 


In  our  assumptions  x  is  large,  while  z  is  finite;  therefore  we  shall 
have  approximately  p  *  -  z. 

In  the  case  under  examination  the  main  part  of  the  integration  in 
integral  (6,02)  will  correspond  to  values  of  t  of  the  order  of  unity.  But  with 
y  large  and  t  finite,  the  asymptotic  expression  (4,07)  will  apply  to  the  func¬ 
tion  w(t  -  y),  which  will  give 


.ii  l  .2,  .NV? 

ir  -  t  My-t) 

w  (t  -  y)  =  e  (y  -  t)  e 


or  approximately 


.  n  1  .2  1  2  .  r 

l4  ~4  "l^t 

w(t  -  y)  =  e  ye 


(6,17) 


(6,18) 


i  u 


Subst  iLuti 
the  multiple 


!  1 

ng  (6,18)  Ln  (6,02)  and  exchanging  the  vaLue  x*y4 
in  front  of  the  integral,  we  shall  have 


for  unity 


v(x»y,q) 


(X  -  /^q>. 


who  re 


l  (  lzt 

Vj(z,q)  =  — —  )  — - - - — —  dt. 

/>  /  w'(t)  -  qw(t) 


(6,19) 


(6,20) 


The  terms  discarded  in  (6,19)  will.,  with  finite  z  values,  be  c.l  the 
order  of  1/r  y  (or  1/x). 

Thus,  in  our  case,  the  function  V(x,y,q)  from  two  arguments  and  from 
Parameter  q  will  come  down  to  the  function  V(z,q)  from  one  argument  and  from 
the  very  same  parameter.  This  is  a  considerable  simplification. 

We  shall  define  the  formula  linking  the  attenuation  function  W  to  the 
value  V.  We  have  the  identity 


2 

3 


% 


(6,21) 


where  m„  has  the  value  (5,18).  Discarding  in  (6,21)  the  last  term,  we 
get  from  (6,03)  and  (6,19) 

i  3 
3Z 

W  =  eJ  (z,q)  (6,22) 


In  this  way,  in  our  approximation,  the  attenuation  function  W  depends 
on  x  and  y  only  owing  to  z  =  z  -  /~y  . 

The  I  unction  Vj(z,q)  is  the  whole  transcendental  function  from  the 
variable  .  For  positive  z  values  the  integral  (6,20)  can  be  computed  as 
the  sum  ol  deductions,  which  gives 
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V]  (•/.,<!)  -  i  2  /v  l  ('t'  -  q'  )w(t's)  (where  55  >  °»  (6>23) 


who re  ts  ~  are  the  previously  investigated  roots  of  the  equation  (6,09). 

Series  (6,23)  converges  more  quickly  the  larger  z  is.  At  large  positive  values 
of  z  its  sum  boils  down  to  the  first  term.  For  finite  negative  z  (e.g.  for 
-  2  ••  v.  <  0)  Integral  (6,20)  must  be  computed  by  squaring.  For  large  negative 
values  of  7.  one  may  compute  the  integral  by  the  fixed  phase  method,  at  which 
one  obtains 


Vi(z,  q)  = 


L  3 

“  3 z 

2e 

1 

z 


(6,24) 


and,  as  :>  result  of  (6,22) 


1  +  — 2- 

z 


(6.25) 


Keeping  in  mind  that  the  approximation  z  =  -  p,  we  get  a  coincidence  with 
Formula  (6,03). 

In  conclusion  we  note  that  our  basic  formula  (6,02)  may  also  be  obtained 
by  a  method  of  parabolic  relation  by  M.A.  Leontovich  and  applied  by  him  to 
the  derivation  of  the  Weyl  -  van  der  Pol  formula.  Use  of  Leontovich's  method 
(with  certain  refinements),  applied  to  this  problem,  will  be  examined  by  us  in 
a  special  article. 


SUMMARY 

The  problem  of  the  propagation  of  radio  waves  around  the  homogeneous 
surface  of  the  earth  is  investigated.  The  diffraction  effects  are  considered 
but  the  ini  1 uencc  of  the  ionosphere  is  neglected.  The  aim  of  the  paper  is 
to  derive  tormulae  for  the  wave  amplitude  as  a  function  of  the  elevation  of 
the  sourie,  its  distance  liom  the  point  of  observation  (situated  on  the  surface 
ol  the  earth),  ol  the  wave  length  and  ot  the  electrical  properties  of  the  soil. 
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The  main  result  is  the  derivation  of  an  expression  for  the  attenuation 
factor  of  an  integral.  This  expression  is  valid  for  all  the  values  of 
parameters,  which  are  of  practical  interest.  In  the  limiting  cases  the  well 
known  formulae  are  obtained:  The  Weyl  -  van  der  Pol  formula  for  the 
illuminated  region  and  the  formula  which  corresponds  to  the  first  term  in 
Watson's  series  for  the  shaded  region  (the  latter  in  a  slightly  corrected 
form).  Essentially  new  is  the  investigation  of  the  region  of  the  continuous 
transition  from  the  illuminated  region  to  the  shaded  one.  Methods  for 
numerical  calculations  of  sums  and  integrals  involved  in  the  problem  are 
elaborated. 
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